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q_j i Abstract. It is a recent observation that entanglement classification for qubits is closely related to 

["t , ' local 5'L(2,C)-invariants including the invariance under qubit permutations [1, 2, 3], which has been 

^_^ , termed SL* invariance. In order to single out the SL* invariants, we analyze the 5L (2, C) -invariants 

^ ■ of four resp. five qubits and decompose them into irreducible modules for the symmetric group ^4 

I resp. ^5 of qubit permutations. A classifying set of measures of genuine multipartite entanglement 

I— I . is given by the ideal of the algebra of 5L* -invariants vanishing on arbitrary product states. We find 

,^ ' that low degree homogeneous components of this ideal can be constructed in full by using the 

Qh| approach introduced in Refs. [2, 4]. Our analysis highlights an intimate connection between this 

.4—) ■ latter procedure and the standard methods to create invariants, such as the H-process [5]. As the 

Ch ' degrees of invariants increase, the alternative method proves to be particularly efficient. 
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*> : INTRODUCTION 

^ \ The quantification and classification of multipartite entanglement takes an important part 

^H ■ in quantum information theory and is subject to a lively discussion in the recent litera- 

t:;^ . ture. Many different sets of measures of entanglement have been proposed for introduc- 

es \ ing some order and insight into the Hilbert space of multipartite systems. An important 

part of the discussion addresses the underlying invariance group the measures have to 
have. Local unitary invariance is clearly a minimal requirement but must be extended to 
. ^ \ invariance under local special linear transformations, when generalized local operations 

/S \ are admitted. For qubit systems this local invariance group is the 5L(2,C). The invari- 

j^ ■ ance group of q qubits is then given by SLq := 5L(2,C)®^, where we will even omit 

the index wherever it doesn't create confusion ' . Interestingly, the demand of invariance 
under SL operations on the measures of pure state entanglement readily implies that the 
induced measure on mixed states is an entanglement monotone when extended through 
its convex roof [6, 7]. The requirement of SL- invariance is restrictive enough to even 
single out a distinguished class of genuine multipartite entangled states: the nonzero 
SLOCC classes are made of all those states that do not vanish after infinitely many SL 
operations (the local filtering operations of e.g. Ref. [7]). Each such nonzero Stochas- 
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^ That is: where either it is clear what number of qubits we are talking about, or in generic statements 
applying to arbitrary number of qubits. We deviate from the standard definition SLq := SL{q,(C); since we 
deal exclusively with qubits throughout the paper, this should not cause any confusion. 



tic Local Operation and Classical Communication (SLOCC) class has a representative 
which can be considered as maximally entangled state within that class. Interestingly, 
for three qubits a maximally three-tangled state with no concurrence exists, and for four 
qubits there are three four-tangled states with neither three-tangle nor concurrence. It 
is an open question whether such representatives carrying only the genuine multipartite 
entanglement classes will exist in general. 

It is worth emphasizing at this point that any function of the pure state coefficients that 
is invariant under SL transformations will remain unchanged by such local filtering oper- 
ations. Consequently, the complementary zero SLOCC class exists that contains all states 
for which all these invariant functions are zero. A prime example for a representative of 
the zero SLOCC class is the multiqubit W state \W) =Zi OCi 10 with |z) = |0 . . . 1 . . . 0) 
being a state with all zeros but a single 1 placed at site number / (or straightforward 
generalizations of it to higher local dimension). Notwithstanding its globally distributed 
entanglement of pairs we therefore would not call it genuinely multipartite entangled. 

5L-invariance has been intensely studied for three qubit systems in Refs. [1,8] and 
for four qubit systems in Refs. [3, 9], and geometric aspects of such invariants have 
been highlighted in Refs. [10, 11, 12]. Preliminary results for five qubits have been pre- 
sented recently [5]. Independent of these approaches, a method based on local 5L(2,C)- 
invariant operators has been suggested with emphasis on permutation invariance of the 
global entanglement measure [2, 4] . Permutation invariance has been highlighted as a 
demand on global entanglement measures already in Ref. [13] and later in Ref. [3], 
where the semidirect product of SLq and the symmetric group Sq of qubit permutations 
has been termed SL* := 5L(2,C)®^ x Sq, which we will abbreviate as SL*. In addition, 
Refs. [2, 4] focus on those invariants that vanish on all product states. These form an 
ideal in the ring of 5L* -invariants and are important for distinguishing genuine multipar- 
tite entangled states from tensor products of entangled states such as \GHZ) ® \GHZ). 
It lies within the nonzero SLOCC class of 6 qubit entangled states but is not genuinely 
6 qubit entangled. 

In this work we will use both local invariant operators as proposed in Refs. [2, 4] 
and the Cay ley Q-process to construct polynomial invariants. First we compare both 
approaches for the known complete set of invariants of four qubits [9] and those for five 
qubits up to degree 6 (see Ref. [5]). For these known cases we follow the notation of 
Refs. [9, 5] and express the invariants presented there in terms of combs and filters from 
Refs. [2, 4]. Then we go considerably ahead up to degree 12 with an outlook to degrees 
14 and 16. 

The manuscript is organized as follows: the next Section reviews the approach to SL- 
invariants using local invariant operators and fixes the notation. Section summarizes 
the main results for 4-qubit invariants from Ref. [9] and establishes the correspondence 
between these invariants and those from Ref. [2]. Section briefly revisits the known 
polynomial invariants (up to degree 6) and then gives a new and complete characteriza- 
tion for the space of five-qubit 5L-invariants for degrees 8, 10 and 12. We also make a 
few remarks regarding the invariants of degree 14 and 16. In Section we determine the 
Hilbert series for 5L* -invariants which confirms the findings of the preceding Sections. 
After presenting an interesting connection between the Cayley ^-process and local in- 
variant operators (combs) in Section we draw our conclusions in Section . The appendix 
provides a detailed discussion of the concepts and notations of the comb based approach 



in Refs. [2, 4] together with a prescription for their evaluation. 

LOCAL INVARIANT OPERATORS AND NOTATION 

Before we start with our analysis, we give a brief summary of the approach using local 
invariant operators. For a more detailed description of this approach see the appendix . 

We will refer to multiple copies of a given quantum state simply as copies in what 
follows. The Hilbert space for m copies of a quantum state x^r of q qubits can be written 
as 

where J^n ® J^2 ® ■ ■ ■ ® J^iq is the Hilbert space for the i-th copy of i//, Mtj is the Hilbert 
space for the j-th qubit in this copy, and • is used as the tensor product sign between the 
Hilbert spaces of different copies of i//. 

We will often use the notion of an expectation value of an operator <9, which for a 
pure state |i//) is defined as (v/| 0|v/). A qubit comb has been defined in Refs. [2, 4] 
as an antilinear operator acting on a single or multiple copy of a pure single qubit state 
{q = 1) which has zero expectation value for all such states. We point out that combs are 
5L(2,C) -invariant operators. Two independent combs 02*^ and o^it^o^ft have been 
identified in terms of the Pauli matrices 

10\ /01\ f -i\ /lO 



^0 •= I 1 ; ' ^1 •= ^ 1 ; ' ""^-[i ; ' ^^3 :^ ^ -1 

' (1) 

where (t is the complex conjugation in the eigenbasis of 03, and the contraction is defined 
via the pseudo-metric G^v '■= ^^vgfi as 
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-- £^AiO-M»f^M (2) 

11=0 

{go,gug2,g3) := (-1,1,0,1) (3) 

Being 5L-invariant, both combs admit the construction of antilinear operators acting on 
multiple copies of pure multiqubit states that are 5L-invariant. Polynomial invariants 
are then constructed from multiqubit operators obtained from combs as their antilinear 
expectation values for a general multiqubit pure state. These invariants are homogeneous 
polynomials in the basis coefficients of the state 1 1//) (see the appendix for more details). 
We will use double brackets to denote the expectation value of an antilinear operator and 
we often omit the tensor product sign (g) : 



((C72C72)) : = 


= {{02 ® 0-2)) := ( V^* 1 (72 <S) 0-2 V^) 


(4) 


o^o\..)) :- 


= (V/*| •(!//* (o-^(g)(7v(S)...)«(a'^OC7^(g). 
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Here, ly/) is a pure state of ^ > 2 qubits (expressible as a vector in (C^)®^). The 
double bracket expressions in Eq. (5) can be evaluated by first calculating the (antilinear) 



expectation values for the two copies of 1 1//) and performing afterwards the contractions 
with the pseudometric. In (4) and (5) we have shown expressions for one and two copies 
of the state only. The extension to more than two copies (and hence higher degree of 
the invariant) is defined analogously. A measure of entanglement is then defined as the 
absolute value of such an invariant, e.g. C= |(((72C72))| is the pure state concurrence [14]. 

By a product state we mean a state that can be written as a tensor product on some 
bipartition of the system. We slightly relax the use of the term "filter" as compared 
with Ref. [2], where also permutation invariance was included. We will call di filter an 
invariant which vanishes on all product states and reserve the term Sh* -filter for those 
which are also invariant under qubit permutations. The algebra of complex holomorphic 
polynomial 5L-invariants resp. 5L* -invariants of q qubits will be denoted by Inv'^^ resp. 
Inv'^^ . The number of qubits, q, will be clear from the context. The filters form an ideal 
J^(f^ of the algebra Inv^^. The 5L* -filters form an ideal J^(f^* := J^q^^ n InV^^* of InV^^* . 
The subspace of Inv"^^, Inv'^^ , etc. consisting of homogeneous invariants of degree d 
will be denoted by adding the degree as an index. E.g. Inv^^ and Inv^^ . 

Let us point out four important facts. First, Inv^^ = whenever d is odd. Second, 
Invj^ = if ^ is odd and has dimension 1 if ^ is even. These facts are special cases of 
more general results proved in [15, Prop. 11.1 and Cor. 11.2]. Third, the dimension of 
Inv^"^ is equal to (2^^^ + ( — 1)^)/3. Fourth, the dimension of Inv4^* is equal to [^-g-J, 
where [t\ denotes the Gauss parenthesis, i.e. the largest integer n <t. These two results 
are proved in [15, Cor. 11.4 and Prop. 11.3]. 

We will also use the notion of a relative SL* -invariant for an 5L-invariant that is 
fixed up to a sign under all qubit permutations. Antisymmetric relative invariants will be 
termed odd relative invariants or 5L!l -invariants. In this context we will sometimes need 
to either symmetrize or antisymmetrize a given invariant for obtaining the corresponding 
SL* and SL^ invariants, respectively. For an operator O, whose dependence on qubit 
permutations is indicated by the permutation operator as an index, we use the definitions 

{0\)a •= 77 E sign;rd;,. 
^- neSq 

An Sq orbit of an invariant / will be denoted by Sq o /, and Kij will denote the permutation 
operator that exchanges qubit numbers i and j. Furthermore we will say that an invariant 
is generically of degree d, if it is not expressible as a polynomial in invariants of lower 
degrees. To simplify the notation, we set V^ ■= Inv^^ and denote by U/j the subspace of 
Vd spanned by the products of homogeneous lower degree invariants, i.e. by VsVds for 
s = l,...,(i— 1. For the decomposition of the space Vj (or Uj) into simple modules Xi 
of the symmetric group Sq we use the notation of Ref. [16]. 



SL AND 5L* -INVARIANTS FOR FOUR QUBITS 

The Hilbert series for SL4 -invariants is [9] 

h{t) ^ 



(l-?2)(l-?4)2(l_^6) 

+ ... 
From theorem 4.2 of Ref. [3] we obtain immediately the Hilbert series for 5L4 -invariants 
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hsL*{t) 



(l-?2)(l-?6)(l_^8)(l_^12) 

= l+?2 + ?V2?^ + 3?^ + 3?^° + 5?^2 + 6?^V7?i^ + 9?^^ 

+ ... 

We deduce that the algebra Inv'^^'* is a polynomial algebra with generators of degree 2, 
4, 4 and 6 and, similarly, that Inv'^^4 is a polynomial algebra with generators of degree 
2, 6, 8 and 12. 

Furtheremore, a complete set of invariants [9] and co variants [17] is known. With the 
focus of finding measures for genuine multipartite entanglement, three independent filter 
invariants have been constructed in Ref. [2] 

^{^ = {{0^0v0202» 0^020x02^020^ 0^02)) 
^2 = (( C^/i Ov 0202^0^ O2 Ox O2 • 

020^020t: • 02020^0^)) 

^3 = - i{0^0v0202 • 0^^0^0202 • Op020T:02 • 

OP 020^02* 0k02020x»0'^ 02020^)) . 

We will use this Section to work out interrelations between the two approaches. 

Degree 2 

As mentioned in the end of Section , this smallest possible degree appears only for 
an even number of qubits q and the corresponding space V2 is one-dimensional with the 
(^-tangle of Wong and Christensen [18] as generator. Here, for q = 4, this generator is 
the 4-tangle and has been termed H in Ref. [9] 

H{y/) ^l{{02020202)) =:^-'rf^''^ . 
It does not vanish on tensor products of 2-qubit entangled states and so it is not a filter. 



Summarizing, we have that 
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Inv2 = span{H} 
0. 



Degree 4 

Besides the one-dimensional space U4 spanned by H^, there exist three new invariants 
of degree 4, namely L, M, and A^ subject to the relation L + M + N = 0. Expressed in 
terms of the coefficients of the wave function 

1 24-1 

|V')= L \ifijki\ijkl)=: Y, an\n), 
with the identification | ijkl) = \i + 2j + 4k + Sl), they are given by the determinants 



ciQ a4 cig a\2 
a\ a^ cig (313 
CI2 ci(f a\Q (314 



M 



qq a% a2 aiQ 

CI4 a\2 Cl(, Cl\4 
as (313 (37 (315 



A^: 



ciQ ai a^ cig 

(32 (33 aiQ Clll 

CI4 (35 (3i2 (?13 

(36 (37 (314 (^IS 



(6) 



They can be expressed in terms of the following invariants obtained from local invariant 
operators 



"^^ 



(4) 



4;(i,2) 

^(4) 
^4; (1,3) 

^4;(1,4) 



{{0^0v0202»0^0''0202)), 
{{0^02(7v02» (7^(72(7" O2)), 
{{On02020v • 0^02020^)). 



Indeed we have 
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and 
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4;(1,2) ^4;(1,3) 
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4;(1,2)^^4;(1,3)^ "4;(1,4) 



For analoguously defined invariants '^^.l ^\, '^4. {2 4] ^^^ '^4(2 3) ^^ ^^^^ ^^ identities 
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4;(1,2) - "4;(3,4) 
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4;(1,3) 



^(4) 
'4;(2,4)' 
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(4) 



^. 
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4;(1,4) - '^4;(2,3)- 



It is interesting to mention at this point that further identities appear besides those 
stated above. Examples are 



{{oyiOvOxO^»o^o''o^o^)) = 36H 



(7) 



and the identity for the three-tangle in [2]. We will report on such identities also for 
five qubit invariants. They suggest that double contractions (o^ (g> (Jv) • {<y^ ® <y^) 
within a pair of copies could be somehow removed. However, the non-trivial example 
{{02020202)) 7^ {{Ojj_Ov0202)) {{0^0^0202)) demonstrates that double contractions can 
not simply be removed. It would be worthwhile analyzing this curious observation in 
more detail and a rigorous reduction scheme would be highly desirable. The interrelation 
between the Q-process and the comb approach in Section singles out one origin for such 
identities. In particular it explains all the identities for degree 4 invariants mentioned 
here. 



Summarizing, we have 



InV4 - Span{<^4.(j2)'^4;(l,3)''^4;(l,4)} 



Inv^^ = span{//^} C U4 



^0-4 = 



Degree 6 

By invoking the Hilbert series, we deduce that dim f/g = 3 . We claim that Vg is spanned 
by U(, and the filter ^| from Ref. [2]. As dimVe = 4, it suffices to observe that 

Defining the 5L*-invariant W := D^y + Dxz + D^t, the expressions for the Duv from 
Ref. [9] give 



H{N-M) -- 


= 3D,y-W, 


H{L-N) = 


= 3Dxz-W, 


H{M-L) = 


= 3D,t-W. 



For comparison with Ref. [3] the correspondence for the invariants are Dja -^ D, D^y -^ 
E, Dxz ^ F and W ^ r. 

The subspace of SL- invariants of degree 6 is spanned by D^y, Dxz, ^xt, and H^ [9] and 
we find that 

^^"^^ = S{4W-H^). 

From these relations all invariants in this subspace are readily expressed in terms of 
comb-based invariants. It is worth noticing that ^j ' is a filter and that it spans the 
subspace J^q.^. 

Summarizing, we have 

Invf = span{//3,//^4(.^)j2)'^^i;a3)'^/'^} 
Invf* = span{//3^^/^^} 
^0^^ = span{^|4)}. 



All spaces Inv^^ with J > 6 are built from generators of degree 2, 4, and 6. What we 
will focus on in the rest of this section is to construct a complete set of generators for 
the ideal ^q^* . 



Degree 8 

In this case we have U^ = V% = Invg^ and dimVg = 7. By proposition 0.1, which is 
proved below, the degree 8 component of ^^^ is only two-dimensional. It is spanned 

by H^[^^ and the symmetrized filter l^'f^ \ . Defining i::=L^+M^ +N^ we find that 

^2^4) = 16(//4 + 4//2(M-L)-16//D;,f-16LM), 

(^f))^ = '-^{SL-H')-'-^H{4W-H^). 

The action of the symmetric group ^4 on the filter ^2 produces three independent filter 
invariants. We have 

Invf = spm{H\H^^^^^2rH^^^^^^y^,H^i'\s40^^'^} 

Invf* = spm{H\H^i'\(^^^'^)} 

.y,f = span{//^(^),(^f)}. 



We note that the orbit 54 o ^2 in the first formula above can be replaced by the three 
invariants: «^1,2))'' Kd,3))' ^"^ ^i;(i,2)<?i,3)- 



Degree 10 

The degree 10 homogeneous component J^(y|o i^ two-dimensional and is spanned by 
H'^^1 ' and H ( ^j ' ) . The last missing ideal generator is obtained from degree 12. 

Degree 12 and beyond 

The 5L* -invariants of degree 12 are to be built from H, W, E and IT := (L - M) (M - 
N){N-L). The filter ^3^^^ = i'^^^.'J^j 2)'^4-'(i 3)'^4-?i 4) ^^ invariant under qubit permuta- 
tions, i.e. it is an 5L*-filter. We find that 

^(4) = -96//2(8E-//4) -64(32n-f-//6) . 



Proposition 0.1 The ideal .^^^ is generated by the invariants AW — H^, 8E — H^, and 
32n + //^ 

Proof: First, it is easy to check that these three invariants belong to ^^^ . Next, let 
/ e J^Q^ be arbitrary. Note that / is a polynomial in the generators H.W.'L and D of the 
algebra Inv^^ . Without any loss of generality we can assume that / is homogeneous of 
degree 2d. The above three invariants can be used to eliminate W, E and D from /. Then 
the corresponding reduced element /q G ^q^ is a homogeneous polynomial in H only. 
Consequently, /o = cW^ for some constant c. In particular, /o must vanish on arbitrary 
product states. Since H however does not vanish on arbitrary product states, this implies 
c = and completes the proof. 



(4) 



Equivalently, the same ideal is generated by the 5L* -filters ^j , ( ^ 
which are functionally independent [2]. This follows immediately from 
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"^f + 18//2^^W^^+48//3^(4)- 



As an important and often cited invariant, we briefly consider the hyperdeterminant, 
Det, of four qubits. It has degree 24 and is given by 



)8o3 



2''3^Det = {H^ - AW)A + (BE - W)B - 4(32n + H 



6\2 



where 



A = 5//'^ + 20iy//^-144E//^ + 16(5^2 -24n)//3 
-96QW1:H^ + 1536E2// + 1921^(3^^ + 8n) , 
B = //^-136E//^-384m/2 + 256E^ 

This can be translated into an expression in terms of H and the filters ^| , ('^2 ) 

and ^3 in a straightforward manner. 

The decomposition of Inv^^ for even J, 2 < J < 12, into irreducible 54-modules is 
given in table 1 . Note that the 5L* Hilbert series confirms the multiplicities of the trivial 
module X\ . 

TABLE 1. The space of polynomial invariants of degree 2 up to 12 into 
irreducible 54-modules. 



degree 




degree 




degree 






2 


Xx 


4 


X1+X3 


6 


2X1 +X3 


8 


3X1+2X3 


10 


3X1+3X3 


12 


5X1+4X3+X5 



It is interesting to briefly focus on specific multipartite entangled four qubit states. 
One prominent class of states is formed by the so called graph states [19, 20]. They 
are created from a fully polarized state in e.g. .)c-direction by successive action of the 
two-qubit entangling operator Uij := ^(1 + O's;/ + (7i,-j — o^-jO^-j) which is also known 
as the control-(73 gate. A complete characterization of graph states for up to seven 
qubits can be found in Ref. [20]. In the case of four qubits, only two graph state 
classes exist. Representatives are the GHZ state and the 4-qubit cluster state [19]. A 
genuinely entangled four qubit state that falls out of this classification, namely \X) := 
^(V2|llll) + 11000) + 10100) + 10010) + 10001)), has been presented in Ref. [2, 4] 

together with an evaluation of the filters <^^ , ^2 ' ^^'^ '^3 on these three states. 
It is worth noticing that the three filters admit for 7 = 2^ — 1 classes of genuine four- 
party entanglement (in the nonzero SLOCC class). Representative states for these seven 
classes of entanglement can be obtained as coherent superpositions of GHZ, cluster, and 
X state from the table of filter values in Ref. [4]. In this sense these three maximally 
entangled states form a basis for the whole nonzero SLOCC class. Consequently, a 
classification of genuine multipartite entanglement in terms of graph states alone is not 
complete. The account of the complementary set of non-graph states as a resource for 
quantum information processing is largely unexplored. 

SL AND 5L-INVARIANTS FOR FIVE QUBITS 

The SL Hilbert series for five qubits has been determined in Ref. [5] as 

^^^ " (l-?4)5(l_^6)(l_^8)5(l_^10)(l_^12)5 (^) 

= l + 5t'^ + t^ + 36t^ + 15t^^ + 22St^^ + 231t^'^ + 1313t^^+1939t^^ 

+ ... (9) 



We have verified the values of the coefficients C2d of t^'^ in Eq. (9) by using the 
formula C2d = (l/(2'^)0E;re52d^(^)^ where x{^) is the character of the irreducible 
representation of S2d corresponding to the partition [d,d] of the integer 2d. This is a 
special case of the formula from Eq. (29) where we replace d with 2d and insert the 
local Hilbert space dimension n = 2 and the number of qubits k = 5. Both numerator 
and denominator in Eq. (8) are even palindromic polynomials of degrees 104 and 136 
respectively. The expanded Hilbert series tells us that there are 5 invariants of degree 
4, a single invariant of degree 6, 36 invariants of degree 8, 15 invariants of degree 10, 
228 invariants of degree 12, etc. In Ref. [5] the invariants up to degree 6 have been 
determined together with 5 invariants of degree 8. 

The first terms of the SL* Hilbert series are (for details see Section ) 

hsj^.{t) = l+t"^ + 4t^ + Ut^^ + 2t^'^ + 39t^^ + 2lt^^ + 130t^^ + n5t^^ + ... . 

In this section we will give a complete characterization of invariants up to degree 12 
and establish a connection with the invariants from Ref. [4]. Since the Hilbert series 
shows that no invariant of degree 2 exists, we start our analysis with degree 4. 



Di: 


= D, 


D2: 


= Dy 


Dy. 


= D, 


D4 


■=Dt 


D5: 


= Du 



Degree 4 

A straightforward calculation shows that the 5 linearly independent invariants Dy of 
degree 4 (v = x,y,zj, u) from Ref. [5] can be written as 

: {{0^02020202 •0^02020202)) (10) 

: {{020^020202 • 020^020202)) (11) 

: {{02020^0202 • 02020^0202)) (12) 

: {{0202020^02 • 0202020^02)) (13) 

^ {{020202020^*020202020^)) . (14) 

Interestingly, these invariants and their generalizations to higher odd number of qubits 
already appeared in Ref. [18]. They form an 55 -orbit, which is nicely seen from their 
explicit forms (10)-(14). The unique 5L*-invariant of degree 4 is then 

It does not vanish on all product states. Therefore, J^q.^ = for five qubits. It is not 
obvious whether this observation can be extended to larger number of qubits. 

The investigation of the full 55 -orbits of a given set of invariants will be a major tool 
for the construction of the complete space of invariants and the determination of the 
5L* -invariants. In the present case, we only needed e.g. Di in order to create all degree 4 
invariants from its orbit. The decomposition into irreducible 55-modules is V4 = Xi +X2. 

It is an interesting consequence of the completeness of (10)-(14) as generators of 
invariants of degree 4 that additional contractions lead to identities. Two examples are 

{{o^OyOx0202 •0^0^0^0202)) = 3{D4+D5)-P, (15) 

{{o^OyO^o.Op.o^'o^'o^o'oP)) = -3P, (16) 

but also note the above mentioned identities for four qubits. Up to the pref actor, equa- 
tion (16) readily follows from the obvious permutation symmetry. 



Degree 6 and singly even q 

The unique invariant F of degree 6 has been created in Ref. [5] by invoking the Q.- 
process. It is an odd function under qubit permutations, corresponding to the irreducible 
55-module V(, = X-j. It cannot be created from the combs. However, see section . Re- 
membering degree 2, where we know from the Hilbert series that Invj^ = for an odd 
numer of qubits, it seems that singly even q is peculiar. This is particularly true as far as 
the comb-based method is concerned. The only invariant of degree 2 that can be created 
from local invariant operators is the ^-tangle [18] {{Oy ^)), which however is identically 
zero when q is odd. This phenomenon draws wider circles, as expressed in 

Theorem 0.1 For an odd number of qubits, all nonzero SL-invariants that can be con- 
structed from combs have doubly even degree. 



Proof: An expectation value {{Oa^Oa2 • • •)) vanishes if it contains an odd number of 
02^ [4]. Since the contraction with the pseudo-metric G^v does not contain 02 this 
implies that for an odd number of qubits there must be an odd number of contractions 
in each copy, leading to an odd number of contractions where all contractions in each 
of the copies are counted separately. Since the copies are always contracted in pairs, 
this number of contractions must be even. Therefore, the number of copies for an odd 
number of qubits has to be even, leading to a degree divisible by 4, hence doubly even. 
This completes the proof. 

It seems that this fact is intimately related to the observed permutation antisymmetry 
of the invariants of singly even degree (we anticipate here that the one-dimensional 55- 
modules of the generic invariants of degree 6 and 10 are spanned by odd functions 
under qubit permutations). The symmetry of the combs under permutations of the 
copies might hinder assymmetry under qubit exchange, even though there is a profound 
difference between copy and qubit exchange. Indeed, as we will see later, nonzero 
antisymmetrizations of comb-based invariants do exist. Nevertheless, it is natural to 
ask for local invariants that are antisymmetric under the permutation of copies; it turns 
out that no such construction exists that connects two or three copies, i.e. there are no 
antisymmetric combs of order two or three. Also notice that no independent symmetric 
combs exist up to degree four ^. 



Degree 8 

We next proceed with a complete discussion of degree 8 invariants. Looking at the 
Hilbert series, the dimension of this space is 36. The 15 products DiDj, I < i < j < 5, 
form a basis of the subspace Ug. This implies the existence of 21 independent invariants 
that are generically of degree 8. Five of these have been constructed in Ref. [5] by using 
the Q-process. We will at first give the decomposition of C/g and Vg into irreducible 
55 -modules and then establish the connection with these five invariants Hy. We find that 

Us = 2Xi-t-2X2+X3, (17) 

Vs = 4X1+3X2 + 3X3-^X5. (18) 

The dimension of Invf^ for 5 qubits is consequently 4, which agrees with the SL* Hilbert 
series. 



The filter 



^f ^ = {{(y^ii 0-^2 ^M3 <^2<72 • a^i a^^ 020^^ 02 • 

0^,020^' 0^^02*0^^02020202)) (19) 



A detailed presentation of this result will appear elsewhere. 



has been introduced in Ref. [4]. We add two new 5L-invariants 

a2aPar(J2aK»a2(J2a^a''a'')), (20) 

^P = 3{{0^0v(7x02(72»(7t:0^0^02(72» 

0^02020pOK»(y^0202aPo'^)) 

+ (21) 

{{a^Ovax02a2 • 0^0^0^0202 • 

Or02020pOK»0^02020Po'^)). 

Notice that the second summand in (21) is an element of Ug. More precisely, 

{{0^0y0x0202^0^'0''0^0202))=3{D4+D5)-P. 

We claim that these two invariants are filters. It is straightforward to check this claim 
for the invariant (20). To see that also (21) is a filter it suffices to show that it vanishes 
on product states. The only partitions that lead to a nonzero value for both terms in 
the above sum are those factoring out either qubits (2,3) or qubits (4,5). The nonzero 
value is a multiple of powers of concurrence and three-tangle (e.g. Cl 2]'^2-(i 4 5)-*' ^^'^ 
the prefactor is found to be independent of which of the two distinct partitions we take. 
It is then straightforward to check that the above combination vanishes also for these 
factorizations, proving the filter property. 

We will show next that, together with P^, the permutation averages of the above 

three filters span the space Inv|^*. The 55-submodule generated by the filter #^| ' has 
dimension 24 and meets t/g in ^n ^2, a 4-dimensional subspace. Thus by selecting 20 
suitable qubit permutations of this filter, we obtain altogether 15 + 20 = 35 linearly 

independent invariants in Vs. To obtain a basis of Vs, we have to add also the filter #g . 
It is interesting that the filter (19) resp. (20) creates a 24-dimensional space of in- 
variants Vg;! rcsp. Vs;2- These two spaces have a 23-dimensional overlap K. Thus 
^8;; = Ti + K; i = 1,2, where 7] are one-dimensional subspaces of Inv^^ . Further- 
more, also the space created from (21), which we will call Vg-^s can be expressed as 
y8;3 = 73 + K", where K C K and the subspace Tj, C Inv'^^ is one-dimensional. Since 
these spaces have been created from filters, the 7] (/ = 1,2,3) are already the elements 
in ^Q^ we have been looking for. These particular invariants are given by 

r,„ := (^P ) ^ , r,„ := (^i" )_ .T,.,o:^ {^f ), 

A detailed analysis of the characters of the resulting irreducible 55-modules leads to the 
decomposition (18). 

We now give the expression of the three invariants in J^q^ in terms of the invariants 
obtained in Ref. [5]. To this end, we define a second 5L*-invariant in Us, namely 

i=l 



and use the sum of all the Hy 



Hn 



LH.. 



i=l 



We find that 



Ti-o = P^-3QeUs, T3;o = Ho + P^- 6Q. 



,SL 



Summarizing, Invg is spanned by P^, T\-fi, r2;o, 73;o, or equivalently by P^, Q, T\-fi, T^-fi. 
The subspace .y^.^ is spanned by T\-fl, T2;o, T^-fi. 



Degree 10 

From the Hilbert series we extract that there are 15 independent invariants of degree 
10, where 5 independent elements of Uio are obtained by multiplying the 5 invariants 

Di, i = 1, . . . , 5 with the invariant F of degree 6. Hence, L'^io = Xg +X7 as an 55-module 

3 

The missing ten invariants are in the 55-orbit of Sf^Q which we construct by invoking 
the following ^-process (for the notation used here see Section ) 



^00222 


= ifj^''' 


^20022 


= ifjr'' 


^20202 


= ifjr''' 


C20222 


= (^20022, -800222) 


^11131 


(y^ r\ 10101 

= l<-20222,/J 


£'20222 


(n /■^01010 
= (^11131, JJ 


-^11311 


(p /-N 10011 
= {E20222J) 


^11111 


= (-Pi 13 11, -820202) 


^10 


= [HuinJ) 



(22) 



^Iq spans a 14-dimensional space which has a 4-dimensional intersection with Uio- In 
terms of irreducible 55-modules the space of degree 10 invariants decomposes as 

yio = X5 + 2X6 + 2X7 



with dimension counting 15 = 5 + 2*4 + 2* 1. This shows that there are no SL*- 
invariants; however, there are two odd symmetric invariants: PF E Uio and the anti- 

symmetrization of Sf jq ' . Both are in the ideal J^q " . 
Summarizing, Inv^^" = J^qio ^ span{P ■ F, (^fg ) }• 



We follow the notation of Ref. [5]. 



Degree 12 

From the Hilbert series we see that the space yi2 of degree 12 invariants has dimension 
228, where a 141 -dimensional subspace Un emerges from lower degrees. The latter 
space decomposes as 

Ui2 = IXi + 10X2 + 8X3 + 5X4 + 4X5 + X6. 

Hence there are 87 invariants that are generically of degree 12. For the complete recon- 
struction and decomposition of this space into irreducible 55-modules we use the filters 
rather than employing the Q-process, since this reduces significantly the computational 
complexity. The origin of this reduction in computational complexity can be understood 
from the analysis in Section . 

We claim that the 55 -orbits of the five invariants 

^n-\ = ((0-Mi<^M2<7M3<720-2«0-^'0-20-2(T;i4a^5« (23) 

020^^020202 • 02020^^0202 • 
0202020^^02 • 020202020^^)) 

^fil = ((<7miO-M2<7m3<^2C5-2«C7'" (7^2(7^4 C72(T2» (24) 

0-/Z5<72a^3c72(7;X6 • 0^^020^^020y,^ • 
<5^,<S2020^,0^ . 0^^02020^0^')) 

^12;4 = ((0-Ml^M2^M3^20-2«C7^'C720-^4 0-^5(72» (25) 

020^^0^^020y^^ • Oy^^020^^020^ • 

O^'^ 02020^^ Oy,^ • 020202020^^)) 

^Sl = ((0-Mi^M2 0'M3 0'20-2«C7'''(J^2a^4 0-2(J2« (26) 

<y^l^O^O^'^0202 • 0^^0^^020^^02 • 
0^0^020^'02 • 0^'>0^'' 0^^0202)) 

'^ill = ((0-Mi^M2 0'M3 0'm4«^M5»«^'''0-2(72CJ2(72» (27) 

020^^020^^0^^ • O^^O^gO^^O^O^' • 

o^^o^'> 020^^02 • 020202020^')) 

generate the full set of 87 new invariants. The symbol ^ indicates that the invariant has 
the filter property. The filters (23) and (24) are taken from [4], whereas the invariants 
(25), (26) and (27) are new. 

We start our construction of a basis of yi2 by choosing a basis of the subspace U12 

(141 elements). Next we make use of the filter ^i2-4- "^^^ ^5 -module that it generates has 
dimension 112 and intersects U12 in a 44-dimensional submodule. Thus we can construct 
the next 68 basis elements by applying suitable qubit permutations to this filter. The next 

15 elements of the basis are obtained similarly from ^12-2' ^^'^ ^ more from ^i2-2- ^^^^ 
gives in total 141+68-t-15-f2 = 226 basis elements. A full basis of Vn is obtained by 
adjoining the invariants ^12- 1 ^^'^ ^12-6 • "^^^^ proves the claim made above. 



It is staightforward to construct filters from ^12-2 ^^'^ "^n-e ^y subtracting suitable 
elements of f/12. In both cases there is a single partition for which the invariant does 
not vanish on corresponding product states: for the partition (1,2) (3, 4, 5) we have 

^12;2 = 9Cf^2 ■ '^3;3,4,5 whereas for (2,4)(1,3,5) we obtain that ^^^.^g = SCfg ■ t|.i 45. 
Filters are constructed by subtracting the [/12-elements 

^2-2 = ((CT^,Cy^2C7^3<72<72«(7^'(7^2(7^3(72(T2))- 

(( (7^5 0-^6 (7^4 O2 02*0^^ (7^7 O2 C7^8 O2 • 

a^,o^^020^'a2 • a^a^' 0^^0202)), 

{{020^^020^^a^,^ • 0^0^020^' 02 • 
0^02020202 • 020202020^^)) 

from ^j2-2 ^^'^ ^i2-6' respectively. 

It is interesting to mention here that an 5L* -filter can readily be constructed from 
combs as follows 

^0^^ = ((0-Mi^Ai2 0'M3 0'M4<^M5»<^^'<^2(720-20-2« 
020^^ O2 O2 O2 • O2 O2 O^^ O2 O2 • 
0202020^^02 • 020202020^^)) 

It turns out that ^q is equivalent modulo U12 to the symmetrization of ^i2i- 

We find that the 7-dimensional space U12 H Invfj is spanned by P^, PTjq (j=l,2,3), 
F^, Y^iD], and /Di^^A . Besides PTj-q (j=l,2,3), (^i^fA and ^^^\ also P^ - 

9E,D3isin^(f^*. 

The complementary 5-dimensional space in Invj2 is spanned by 



'^12;!/,' \^12;2/.' \^ V1;A j ' \=^12;2/^, ' \=^12;6y ^, 



r(5)\ /^(5)\ /^(5)\ /^(5) 

12;2 

The two antisymmetrized filters 



^\l\ll ^ ' \^12;4 



span the space of 5L^ -invariants of degree 12, whereas / 3^j2-6 ~ ^12-2 ) ^^ ^'^ '^i^ ■ ^'- 
is worthwhile noticing that the comb-based invariants also create the SL!L -invariants in 
V\2\ those of degree 6 and 10, which are not accessible by the comb approach, are not 
needed. 




Summarizing we have a 10-dimensional space ^^.^2 inside a 12-dimensional space 
Invj2 • In addition there are two 5L!l -invariants, both belong to the filter ideal. Thus, 

The complete decomposition into irreducible 55 -modules is given in table 2. 

TABLE 2. The space of polynomial invariants of degree 2 up to 12 into 
irreducible ^5 -modules. 

I degree I 11 degree I I 

I 2 I II 4 I Xi +X2 I 

I 6 I X7 II 8 I 4X1 +3X2 + 3X3+ X5 I 

I 10 I Xs+2X(, + 2X-i II 12 I 12X1+15x2+14x3+6X4+8X5+2x6+2x7 I 



Beyond degree 12 

We add here a couple of remarks about the next two degrees, 14 and 16. To this 
end, let Y denote a minimal set of (homogeneous) generators of the algebra Inv^^. We 
know that 7 is a finite set, but its cardinality is not known. It is a disjoint union of 
the subsets 7^/ := 7 fl Inv^^. From the Hilbert series we know that |y^/| =0 for odd 
d and for d = 2. Our computations show that for d = 4,6,8, 10, 12, 14, 16 we have 
17^1=5,1,21,10,87,145,247. 

Let us assume that the conjecture made in Ref. [5] regarding the Cohen-Macaulay 
ring structure of Inv'^^ is correct, i.e., that the primary invariants consist of five polyno- 
mials of degree 4, one of degree 6, five of degree 8, one of degree 10, and five of degree 
12. Then Inv"^^ would be a free module of rank 3 014 400 over the algebra generated 
by the primary invariants (a polynomial algebra in 17 variables). Moreover, the coeffi- 
cients of the numerator of the Hilbert series give, for each degree, the number of basis 
elements of this free module. The first six nonzero coefficients are 1, 16, 9, 82, 145, 383 
and the degrees of the corresponding basis elements are 0, 8, 10, 12, 14, 16, respectively 
(see Table 1 in Ref. [5]). For instance, for J = 16 we have 383 basis elements. We may 
assume that this basis contains 7i6. Consequently there must be 383 — 247 = 136 basis 
elements of degree 16 that come from the products of basis elements of degree 8. As 
there are sixteen basis elements of degree 8, the number of different products of two 
of them (including the squares) is indeed 136. This may be interpreted as additional 
evidence for the validity of the above mentioned conjecture. 

It is interesting to examine again the class of graph states. For five qubits there exist 
four inequivalent graph states: numbers 5 through 8 in Ref. [20] (see figure 1). They are 
distinguished already by the degree 4 invariants D,, / = 1, ... 5. Among 5L* -invariants 
up to degree 8, the state d) of figure 1 is only detected by Tqjs (and hence also by Hq). It 



a) b) c) d) 







FIGURE 1. Full set of inequivalent graph states on 5 qubits as graphically depicted in Ref. [20]. a) is 
equivalent to the GHZ state |0) + |31), where the binary decomposition of the numbers gives the qubit 
state, e.g. |31) = |lllll).b)is equivalent to |1) + |2) + |28) + |31) and emerges from the 4-qubit cluster 
state after copying and stacking one qubit (telescoping), c) is equivalent to |1) + |6) + |24) + |31) and 
emerges from the 3-qubit GHZ state |1) + |2) + |4) + |7) after two-qubit telescoping. Finally, d) is the 
superposition of two orthogonal and telescoped 4-qubit cluster states. They are distinguished already 
by the invariants D, of degree 4. However, all D, vanish for state d). The 5-tuple of (relative) SL* 
generators up to degree 8 - (P,F,Tio,T2o,Tj,q) - evaluated on these four states gives a) (5,0,-1,1,4), 
b) (2,0,^,-i,-|),c) (-1,0, jL, ^1,8), andd) (0,0,0,0,6). 



is worth mentioning that the 5L!l -invariant F does not detect any of these states. Besides 
the two maximally entangled states from Ref. [4] 1 1) + |2) + |4) + |24) + V2|31) and 
W4 + \/3|31), also superpositions of these two states and possibly including the four 
graph states fall out of the graph state classification. This basis of 6 states would thereby 
admit for already 2^ — 1 = 63 distinct SLOCC classes of 5-qubit states. Only four of 
these are covered by graph states. 

CHARACTER COMPUTATIONS AND THE HILBERT SERIES 

FOR SL* 

It is interesting to mention that the coefficients of the Hilbert series for the symmetry 
group SL and SL* can be obtained directly using the results of Ref. [15]. Here, we recall 
some results from that work and use them to compute the dimension of the space of 
5L* -invariants of degree 2d. We also do the same for the SL^ invariants. 

As in the cited reference, we shall be more general. First, instead of qubits we may 
work with qudits, i.e., we consider the vector representation of SU{n) or 5L(n,C) on 
V = C". By taking k copies of 5L(n,C) and k copies of V and tensoring, we obtain the 
standard representation of G = SL(n, C)®^ on V®^. Let us denote by Mn,k the algebra of 
holomorphic polynomial functions on V®^, and by ^n,k.d its subspace consisting of the 
homogeneous polynomials of degree d. Next, denote by ^^^- ^i the subspace of ^n,k,d 
consisting of G-invariant polynomials. If d is not divisible by n, then ^^^ ^ = by [15, 
Proposition 11.1]. 

Assume now that d = nr and let K = [r"] be the partition of d into n equal parts. 
Denote by Ej,; the irreducible module of the symmetric group S^ which corresponds to 
n, and let x be its character. Then by [15, Corollary 11.1] we have an isomorphism 

^^Xd=(EfY'' (28) 



of 5yt-modules. The superscript S^ means that one has to form the space of invariants of 
Sj, i.e., the largest trivial S^ submodule ofE®'^. The other symmetric group, Sk, acts on 
both sides by permuting the tensor factors. 

This formula is very useful. For instance, one obtains immediately the following 
formula for the dimension of the space of G invariants of degree d: 

dim<,,, = jy£z(^)'. (29) 

By symmetrization, i.e., by taking the 5,t-invariants on both sides of Eq. (28), and 
taking into account that the actions of Sk and S^ commute, we obtain that 

as complex vector spaces. (By S^{Ejt) we denote the k-th degree piece of the symmetric 
algebra S{Ejt) of the module E^.) 

By performing anti-symmetrization instead of the symmetrization, one obtains a 
similar formula for the dimension of the space of odd invariants of Sk in ^^kd- Then 

on the right hand side one should replace the symmetric power S^{Et^) by the exterior 
power a'^{Eji). 

The character x^ '^ of the 5^/-module S^{Eji) is given by the classical formula [21, 22] 



k 



(i) a=l '«-^ 



(30) 



where the summation is over all sequences (i) = {iij2i • • • , h) of nonnegative integers 
such that 

J^ a /« = k. 



a 



This is valid for all permutations g G Sj. 

Similarly, the 5^-character x of the k-th exterior power A^{Eji) is given by the 
formula 






#'M-En ^"'.-.jr ' ■ (31) 

(i) a=l 



The values of the irreducible characters of Sd are easily available, say in James and 
Kerber book [16] or in software systems such as Maple or GAP. Hence, we obtain the 
following formula for the space of joint G and S,t-invariants of degree d = nr. 



dim(<,,)'^ = ]^Iz(^)(,). 



In our case we have n = 2, since we work with qubits, and k = 5, i.e., the number 
of qubits is 5. In that case there are exactly seven sequences (i) having the required 



property. Explicitly, they are: (5,0,0,0,0), (3,1,0,0,0), (1,2,0,0,0), (2,0,1,0,0), 
(0, 1 , 1 , 0, 0) , ( 1 , 0, 0, 1 , 0) and (0, 0, 0, 0, 1 ) . Formula (30) then reads as 

+20x{g^)Xig')+^0x{g)xig^) + 2Axig'). 

For instance, if J = 8 = 2 ■ 4 we have r = 4, ;r = [4, 4], Ej^ is the module Xg, and the 
values of x on the representatives of the 22 conjugacy classes of Sg are 

14,4,2,0,6,-1,1,-1,2,-2,-2,0,2,1,2,-1,-1,-1,0,0,0,0 

(see James and Kerber, p. 351). By using the above formula, we find that the values of 
the character x^^' on the same representatives are 

8568,216,72,0,536,0,0,0,18,0,-12,0,12,0,24,3,1,0,0,2,0,0. 



,v^,v^,vy, -i-v^,v^, j-^,v^, j-^,v^,^ • T'J ^7'-'7'-'7-^7"7' 



Then the multiplicity of the principal character (i.e., the character of the trivial Ss- 
module Xi) in x^^' is equal to the dimension of the space of 5L* -invariants of degree 
8. Hence we have 

dimlnvf* = - Y^X^^Hg)- 

The evaluation of this sum confirms our finding that this dimension is 4. 

In conclusion, we summarize the results of our computations. The number of linearly 
independent 5L* -invariants in degrees 0, 2, 4,..., 22 is 1,0, 1,0, 4, 0, 12, 2, 39, 21, 130, 
115 respectively. The number of linearly independent relative 5L* -invariants in degrees 
0,2, 4,..., 22 is 1,0, 1, 1,4,2, 14, 11,49,58, 185, 269 respectively. 

CONNECTION BETWEEN THE ^-PROCESS AND THE COMB 

APPROACH 

In this section we present a rephrasing of elements of Cayley's Q-process in terms of 
local invariant antilinear operators. The central operations in the Q-process are determi- 
nants of derivatives 



Q.X = det 



(9/ (9/ 
djf dji 

Xq Xj 



(32) 



with subsequent "trace" tr : x ,x -^ x applied to functions of the wave function coef- 
ficients \]/i^„,,^i^ dressed with auxiliary variables z] such that a wave function |^) := 

L ¥ii,-,iq I '1 , • • ■ , iq) is mapped to the function / := L Viu-,iq4i ' ' ' % ■ ^ typical step 
in the ^-process is then prescribed as [5] 

(p,e)^-'^-=tra%...aV(zOe(/) . 

The key observation is that the action of Q.x (32) amounts to a contraction of two of the 
wave function coefficients with the antisymmetric tensor Ej^i, kj E {0,1} with £oi := 1 



'^. We illustrate this procedure in the most simple example 

-022200 — [J, J) 

— H'li,...,l2,k,lYji,...,j3 Z,-j Z,2 ^13 ^jl ^J2 ^h ' 

5222ool2(.)=i = ¥iu-,i3,k,i¥h,-,h^'' = -((3Jaa2C72)) , 

where we used Einstein sum convention and contraction via e. The above so-called 
transvectant 622200^ which is bilinear in the z^-'' (J = 1,2,3), coincides with the subse- 
quently shown antilinear expectation value after setting all z} = 1 ; here J = I , , 

It is seen that the pairs of wave function copies to be contracted with £ by Q^ (32) are 
all those, whose functions / contain the variable z' and z", respectively. Since both these 
variables might occur in more than one function /, the action of Q.^ will in general lead to 
a sum over such e contractions involving different pairs of copies of the wave function. 
Though it is clear now that each invariant constructed by the ^-process can be directly 
transcribed into a sum of complete contractions of the wave function coefficients via 
the antisymmetric tensor £ = 102, it cannot be directly written in terms of antilinear 
expectation values of 02. A simple three qubit counterexample is the invariant 

whose absolute value is the three-tangle [13]. For obtaining the second equality in 
Eq. (33) we grouped in pairs the first and the last two wave function coefficients. The 
contractions of the idexes ai, a2, b\ and b2 then appear inside these pairs and we will 
call them inner contractions. On the other hand the contractions of the indexes a^ and 
a^ involve two different pairs of coefficients, and we will call them cross-contractions. 
As stated above, for more complicated invariants produced by the Q-process, this 
correspondence is not given by a single full contraction; however, each of those complete 
contractions contained in the Q-process is an invariant. E.g. the invariant F of degree 6 
(see Ref. [5] for its construction from an Q-process) is equivalently expressed as 

Wi . . i&,hi -^fkjhJu . . yifi9,'i2,ho,'i4,h5 

T ^25^3 r ^14)^15 T 

Since the space of degree 6 invariants for five qubits is one-dimensional, the expression 
(34) reproduces precisely this unique invariant F (up to a prefactor). For the invariant of 
degree 10 as constructed in Eq. (22), a possible transcription is 

C#(^) — i;/ 11/ i;/. h,hM. yj6,h . . '5 (!.';.) 

S'lO — Ylul2,l3,H,l5 Yl6,llM,'9,nQ Ym '12 Y ll},lU V-JJ^ 

.i/'l . '8. '10 .,/•. '15. k. iir'n. . '14, '12 
Y '15 '16 Yin '18 '19 r '20, '21 



.//■'22,'20,'13,'16,'25 ii/-. . '18 . . ii/-'17,'23,'21,'24,'19 
Y r'22,'23 '24,'25 Y 



^ It is worth emphasizing that £ = iOi has the physical interpretation of a spinor-metric [23]. 



The space of degree 10 invariants has dimension 15, and the tilde indicates that the 

expression (35) cannot be expected to coincide with ^jq as created from the Q-process. 
The latter is rather a sum over all possible £-contractions emerging from the given Q.- 

process, and ^j^^q is only one element of this sum. This transcription therefore already 
reduces the computational complexity of such invariants. Interestingly, the symmetric 

group 55 generates from $fjQ a 14-dimensional subspace where only PF is missing to 

give the whole 15-dimensional space Vjq . 

A view back onto Eq. (33) suggests a connection between the Q-process and the 
invariant construction from combs; namely that the cross-contraction on the third qubit 
might be substituted by the comb of second order a^ • a^, possibly including a term 
proportional to 02 • O2. In order to make this connection a rigorous statement we 
translate the index contraction into an antilinear expectation value. The symmetry of 
antilinear expectation values 

(VA*|A|(p) = ((p*|At|v/) = ((p*|A|v/) (36) 

for Hermitean operators is crucial for this to work. The procedure is best explained 
graphically in figure 2. The three-qubit wavefunction coefficients are sketched as a staple 





c^^ P c^^ c^ 



FIGURE 2. Three-qubit wavefunction coefficients sketched as a staple of three circles. Each contraction 
with the antisymmetric tensor e = 102 is visualized by a line connecting two circles. Arranged in pairs, an 
expectation value with e corresponds to an inner contraction. 

of three circles there. They are arranged in pairs, and an expectation value corresponds 
to an inner contraction of that pair - i.e. a contraction of wavefunction indexes inside 
such a pair. Each contraction with the antisymmetric tensor e = i02 is visualized by 
a line connecting two circles. The cross-contractions are not yet expressed in terms 
of expectation values. Fortunately, suitable permutations of copies, which are however 
local in the qubits, do exist as to transform also the cross-contractions into expectation 
values without disturbing the inner contractions. In the following we describe this 
iterative procedure. The first equality in Fig. 2 is due to the symmetry (36) for antilinear 
expectation values of Hermitean operators. The second equality is formally expressed as 

3 3 

(V/*|«(V/*|P(72(72(72»0-2(72(72|y)»|V^) = (V^* | • (V^* | (72*72(72 • 0-2a2a2P 1 1//) • 1 1//) 

where P is the symbol for a copy permutation operator and the number three on top of P 
indicates that this permutation operator acts non-trivially only on the third qubit. Using 
P = ^ Lu=o <7^ • <7;i, a straightforward calculation produces 

(CJ2«C72)P = M^v(7^»(7v ,^7) 

= -5 ((7^ •(7^ -(72 •(72) 





where M^v = 5^v"^/i, {'nQ,mi,m2,m^) = (1,— 1, 1, — 1)/2. The resuhing antilinear ex- 
pectation value of Mp^y Op, • Oy is then indicated graphically by a double line connecting 
the copies. For completeness we mention that 

(a^ • a'^) P = -- (c7^ • C7^ + 3a2 • (72) 

which readily follows from the identity {02 • 02) PP = {02 • 02) and Eq. (37). 

In order to see in how far every invariant constructed with an Q-process can be ex- 
pressed in terms of expectation values of local invariant operators, the identity described 
in Fig. 3 is helpful. 



C^^) (^^) ©=^) (^^) 

o — e C3 — e 

FIGURE 3. Using copy permutation operators that act separately on demand on specific qubits, the 
above graphical expression for an identity is obtained. It shows that the rephrasing in terms of antilinear 
expectation values can indeed be obtained in this way. 

It means that the translation of the ^-process into antilinear expectation values can 
be performed iteratively qubit per qubit: an apparent incompatibility of a contraction 
pattern with a fixed ordering of wave function coefficients (vertical rows of circles) is 
resolved iteratively making use of the symmetry (36). It is worth mentioning at this point 
that grouping the wavefunction coefficients in pairs is arbitrary and therefore we obtain 
the same invariant when changing this order by permuting the order of the wave function 
coefficients. This however changes the contraction scheme for the invariant and leads to 
identities for invariants as observed above. For degree 4 invariants this freedom means 
that interchanging 

02*02^ -{02 •02-ap» a^) 

on each qubit leads to the same invariant. This immediately leads to the identities (7), 
those mentioned just above it, and the identities (15) and (16). 

For higher degree invariants, more such pair permutations may occur. Since the 
symmetric group is generated from nearest neighbor exchanges TTjj+i and by virtue 
of the relations 7Ci2^23^n = ^23^12^23, ^h = X it is sufficient to consider the results for 
up to three permutation operators. We find that 

(a2»C72»0'2)Pl2P23 = I [<72 • (72 • (72- 

[Op • o^ • 02 + Op»02» 0^^ + 02* Op • a^) (38) 

—i£klmT^k*T^l*T^m\ 



and 



((72 • (72 • O2) P23P12 = i [(72 • (72 • (72- 

[Op • a^ • a2 -H a^ • (72 • (7^ + (72 • (7^ • (7^) (39) 



where Ti := ao, T2 := Ci, T3 := (73. A further permutation acting on the antisymmetric 
part iEMm'^k*^ • '^m Icads to terms of the type (c^ • c'^ • (72 + cr^ • (72 • (7^ + (72 • O'^ • O''^) 
and vice versa. 

Besides an antisymmetric part in the exchange of copies, only O2 and o^ • a^ appear 
and the set of locally invariant operators is closed. The antisymmetric part is not cap- 
tured by the two combs but in our analysis it appeared to be irrelevant for the search of 
5L* -invariants. It leads, however, to invariants that are antisymmetric under qubit permu- 
tations (see e.g. the invariants of degree 6 and 10 for five qubits). Since an entanglement 
measure is defined as the modulus of an invariant, the requirement for a class-specific 
global entanglement measure must be relaxed to include also 5L!l -invariants. 

CONCLUSIONS 

We have presented a thorough analysis of the polynomial SL-invariants of four and 
five qubits with particular emphasis on the filter ideal, i.e. on those invariants that 
vanish on all product states. Within the complete ring of invariants the filter ideal 
plays an outstanding role because only these invariants can clearly discriminate the 
genuine multipartite entanglement from that within parts of the system. It therefore hosts 
candidates for a class specific quantification of genuine multipartite entanglement and 
its knowledge plays a key role for a systematic analysis and a deeper understanding of 
the structure of entanglement in multipartite systems. A reasonable measure for global 
entanglement should also be invariant under qubit permutations. In order to determine 
the dimension of the subspace of permutation invariant elements in the filter ideal, 
we analyzed the decomposition of the space of polynomial invariants into irreducible 
modules of the symmetric group. 

In the case of four qubits, the standard approach from invariant theory, employing the 
well established Q-process by Cayley, has already lead to the construction of a complete 
set of 5L-invariants [5] . We have compared this approach to an alternative proposal based 
on local invariant operators, termed combs [2, 4]. We could demonstrate that also the 
latter approach generates a complete set of invariants, and we provide a full dictionary 
for expressions from both approaches. We have furthermore established the equivalence 
of the Q-process and the contraction with the spinor metric and provide the missing 
link between the Q-process and the construction from local 5L- invariant operators. This 
implies that major part of the ring of invariants can be generated using the comb based 
method. Indeed we find that all 5L* -invariants and even many 5L!l -invariants in this 
work are accessible to the comb based approach. In addition, the interrelation between 
^-process and the comb approach implies the existence of interesting identities among 
sets of invariants and readily explains those identities observed among invariants of 
degree 4. 

We single out two major advantages of the approach from local invariant operators. 

I Control: The comb based approach admits a high degree of control over specific 
important properties of the invariants that are to be constructed. Of particular 
relevance is the ab initio knowledge about the set of product states for which the 



invariant will vanish. This is a key quality that admits a systematic construction of 
filter invariants; this provides a targeted construction of the filter ideal. In contrast, 
from the Q-process, and equivalently from the contraction with the spinor metric, 
there is no a priori knowledge about the invariant's value on product states. 
II Computational complexity: The interrelation between Q-process, £ contractions 
and the comb based approach explains the notable difference in computational 
complexity we observed. It is clear from these interrelations that the computational 
complexity of the comb based approach is significantly lower than that using 
contractions with the spinor metric e, which itself already constitutes a speed-up as 
compared to the Q-process. This discrepancy grows more important with increasing 
degree of the invariants. 

These advantages permit us to go considerably further in a thorough analysis of in- 
variants, and we demonstrate this for five qubits: we give a complete analysis of SL- 
invariants up to degree 12, and provide an outlook on the situation for degrees 14 and 
16. Although the five qubit case is still far from being completed, we have presented 
a straightforward technique for how to proceed; we are confident that a minimal set of 
generators can be obtained in the way described in this manuscript. All results have been 
cross-checked with predictions from the Hilbert series. To this end we also derived the 
first terms of the Hilbert series for relative 5L* -invariants. 

We hope that the high degree of control paired with the significantly lower compu- 
tational complexity in the generation of invariants have future impact in both quantum 
information and invariant theory. Further analysis would be necessary in order to find an 
expression of all 5L!l -invariants in terms of (antilinear) expectation values. It would be 
also interesting to extend an analysis along the lines proposed e.g. in Refs. [12, 11] in 
order to see whether the filter ideal has a distinguished geometrical interpretation. 

COMB-BASED INVARIANTS 

In this appendix we give a detailed elucidation how comb-based invariants are calcu- 
lated. 

Let the pure ^-qubit quantum state |i//) be expressed in terms of a basis ^ made 
of tensor products of eigenstates | — 1) and |1) of the Pauli spin operator o^, such that 
Ot, \s) = s\s) for 5 = ±1. That is, we have 

^= {\si)^---^\sg) : Sj = ±l} 

In this basis the Pauli spin operators (consider ^ = 1 for the sake of simplicity) assume 
the matrix representations o^''^ := (^'l Ot \s) as given in Eq. (1). Matrix elements of q- 
qubit operators are then defined in the standard way for arbitrary (/-qubit pure states | cp), 
\\f/)as 

((p|ai,®---®oi JV/) := {cp,.0*<yi:'''---<y-r'>su...,s, (40) 



Vj,...,i-^/ n iq 



\(Ecp\,,a:-''...a:fW.u-,s, (41) 



within Einstein summation convention, and [CC<p]^i 



V5'*^ 



(■^1, 



■T^q 



Cicp). 



In this sense the antilinear expectation values as defined in Eq. (5) are specific matrix 
elements of an antilinear operator A = L^^- Here, C is the complex conjugation as 
defined in Eq. (41), and La is the linear operator associated to A. In all this work, the 
operators A and La are antilinear Hermitean and Hermitean, respectively. In the case of 
a single copy of the state we then have 

{{La)) := {W\A^WY = {w\A\wT = {W*\La\w) 

which is a matrix element as defined in Eq.(40) where |(p) -^ |i//*) (see Eq. (41)). For 
any indexes /i , /2, • • • , '^ £ {0, 1 , 2, 3 } we therefore define a bilinear form 

whose value at ((p, i//) is the multiple sum (using the Einstein convention) 



ai,bi a2,b2 



^^a^bn 



(Pai,a2,-,aqWbub2,-,bg- 



This can be also expressed as 



{<^h<^i2 



■(yi,){(p,¥) = {(p*\^ii^(yi: 



'(^M- 



As the first example, we set ^ = 1 and ?i = 2 and we obtain the 5L(2,C) -invariant 
bilinear form 



(<^2)(<P,V^) 



V^i ¥2 



(42) 



However, in this case we have {{02)) '■= ((72) (v^, i//) = for all 1^, which is the comb 
property of the operator 02. 

As another example we take q = 2 and /i = (2 = 2. Since 02 = —ie, we have 

((72a2)((p,V/) = -£""'' £''''''(Pa„a2¥bub2 



<P2,1 <?>2,2 
¥l,i V^l,2 



<Pl,l <Pl,2 
V^2,l V^2,2 



an SL-invariant bilinear form. In the case when 9 = 1//, we obtain the nonzero SL- 
invariant quadratic form 



((a2C72)) := {0202){¥,¥) 



¥i,i V^i,2 

¥2,1 V^2,2 



For operators acting on m copies of the state just replace |i//) by |i//)«-- ••|i//) =: Iv^)"" 
and the corresponding expression for (i//* |. To outline this in more detail, let J^ denote 
the Hilbert space of a single qubit, and J^q — J^®^ the one for the system of q qubits. 
We shall also use the Hilbert space for m copies of this multipartite system. In that case 
we use • to denote tensor products of Hilbert spaces of different copies. Let us now take 
a collection of m bilinear forms of the above type. 



l^^iW^/jW 



■^/,w)' 



1,2,, 



,m. 



and let us form their product 

which is a bilinear form 

The value of this bilinear form on the special elements 

is equal to 

m 

If we further specialize (p^^^ = i//^^^ = i// for all k, we obtain the 2^-ary form of degree 
2m in the complex components of i//: 

•m m 

((11 ^.i« • ■ ■ ^L«)) •= n(^/i« ■ ■ ■ ^lw) (v^' v^) • 

k=l k=l 

We refer to this form of degree 2m as the associated form of the bilinear form (43). This 
definition extends immediately to any bilinear form (44). 

In general, the forms of degree 2m constructed above are not 5L- invariant, but we can 
use their suitable linear combinations to obtain SL- invariant forms. 

In order to do that we proceed as follows. First we select a site, say 5, 1 < 5 < ^, of 
our multipartite system and choose two different copies of the state, say copies p and q, 
I < p < q < m. Next we replace in (43) the Pauli matrices o (p) and a. (,) with symbols 
o^ and cr^, respectively. This is to indicate that the two idexes /i are to be contracted by 
using the pseudo-metric G^y We now interrupt our description to give an example. 

When q = 3 and m = 2 the expression (43) has the form 

We now choose p = I, q = 2 and s = 1. By applying the above instruction, we obtain 
the expression 

(a^a(i)a(i))(c7'' 0.^(2)^(2)). 

By fixing o (i) = cr. (i) = o (2) = c (2) = 02 and performing the /i-contraction this gives 

l2 13 '2 ^3 

the linear combination 

-(o-oO-2(72) (O'0C72O-2) + (<7i(72<72) ((7i(72<72) + ((730-20'2) ((73a20-2)- 

The associated quartic form is then obtained as 

3 

{{0^0202 • 0^0202)) •■= Y. g^i {¥* I o-^fy2(y2 1 vf 

11=0 

3 

V^.'i ,4,4 V ^2 ^2 V^^l,^2,^3 






It generates the 5L* -invariants for three qubits; its modulus is the three-tangle [13]. 

To continue our description, we choose a collection of triples (5,, /?,, ^,), / = 1 , 2, . . . , ? 
such that I < Pi < qi < m and whenever si = Sj, with / 7^ j, we require that the four 
integers p,, qt, pj, qj be all distinct. For each index i, we replace the Pauli matrices on the 
site Si and copies p, and qi with the symbols o^,. and a^', respectively. Next we replace 
all other Pauli matrices in (43) with the matrix 02. Finally, by using the pseudo-metric 
Gj^v, we perform the /i; contractions for each z, I <i<t. We obtain a linear combination 
of bilinear forms of the type given by (43). We refer to these linear combinations as 
comb-based bilinear forms. 

These comb-based forms are homogeneous multilinear expressions in the (complex) 
state coefficients, which are SL^-invariant. This invariance harkens back to the 5L(2,C) 
invariance of the antilinear single qubit combs 02^ and c^C^* o^<t. We formulate this 
statement in 

Theorem .2 Any comb-based bilinear form (and, consequently, also its associated form) 
is an SL-invariant. 

It has been stated in [2, 4] that the combs are SL-invariant, but there is only implicit 
reference to the fact that this derives from the central comb property to have zero 
expectation value on all the local Hilbert spaces. Here we sketch a proof for this 
connection. 

Proof: The comb property for the operator 02, namely that {o2){^f^^f) =0 for all 
single qubit states i//, can be read off directly from Eq. (42), and it can be checked by 
direct calculation that it is the unique operator with this property up to rescaling. Also by 
direct calculation we find that (c^ • c^ )(!//• i//, !//• i//) = for all single qubit states i//. 
Furthermore, this is the unique form (up to rescaling) on M'*"^ satisfying this condition 
which is symmetric under copy-permutation and orthogonal to {02 • 02) in the sense of 
the vanishing scalar product tr {02 • 02) ■ {o^ • o^) = 0. For arbitrary S G 5L(2,C) we 
then find that 

0={O2){SXj/,SXI/) = {S'02S){XI/,XI/) 

for all i/A and, due to the uniqueness property for the operator 02, this implies S'o2S = 02. 
Analogously we have 

= (a;x«(7^)(Sv/«5vA,5vA«Sv/) 

This proves that the two comb operators are 5L(2,C) -invariant. Consequently, a (^r-qubit 
form constructed from those is seen to be SL^-invariant by wrapping a transformation 

S*-^-* —Si®---®Sq with Sj G 5L(2,C) back onto the states. This completes the proof. 

We refer to the 5L-invariants constructed in this manner as the comb-based invariants. 
In some cases this invariant may be zero. 
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